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ABSTRACT 

The filament IRDC G035. 39-00.33 in the W48 molecular complex is one of the darkest infrared clouds observed by Spitzer. It has 
been observed by the PACS (70 and 160 jum) and SPIRE (250, 350, and 500 yt/m) cameras of the Herschel Space Observatory as part of 
the W48 molecular cloud complex in the framework of the HOBYS key programme. The observations reveal a sample of 28 compact 
sources (deconvolved FWHM sizes <0.3 pc) complete down to ~5 M© in G035. 39-00.33 and its surroundings. Among them, 13 
compact sources are massive dense cores with masses >20 M©. The cloud characteristics we derive from the analysis of their spectral 
energy distributions are masses of 20 - 50 M©, sizes of 0.1-0.2 pc, and average densities of 2 - 20 x 10^ cm~^, which make these 
massive dense cores excellent candidates to form intermediate- to high-mass stars. Most of the massive dense cores are located inside 
the G035. 39-00.33 ridge and host IR-quiet high-mass protostars. The large number of protostars found in this filament suggests that 
we are witnessing a mini-burst of star formation with an efficiency of ~15% and a rate density of ~40 M© yr~^ kpc~^ within ~8 pc^, a 
large area covering the full ridge. Part of the extended SiO emission observed towards G035. 39-00.33 is not associated with obvious 
protostars and may originate from low-velocity shocks within converging flows, as advocated by previous studies. 

Key words. ISM: cloud - Submillimeter: ISM - Stars: formation - Stars: protostars - ISM: G035.39-00.33 - ISM: W48 



1, Introduction is the W48 molecular complex. This survey will ultimately 

^ ^, T , . jttt; ri u^.^ i . ^ provide an unbiased view of the formation of OB -type stars 

C3 ^ The Herschel observatory ^Pilbratt et aL||2010D, operatmg at tar- influence of the ambient environment on that process 

infrared to submilimeter wavelengths, is an excellent tool for ^^^tte et al. 2010; Schneider et al. 2010b; Hen nemann et alJ 

studying the early phases of star formation and the connection |^ |di Francesco et all i2010t iHilletaU |201 1). A picture is 

of star precursors to the ambient cloud To address fundamen- ^^^^^ ^ ^^^^ high-mass stars are formed from more 

tal questions of massive star formation, the guaranteed-time key dynamical processes than low-mass stars. Massive dense cores 

P'°^f™t,.fcra;'u 'r^'"^ survey of OB Young Stellar ob- ^^^^s, -0.1 pc, cm'^) either IR-bright or IR-quiet, where 

jects (H0BY5O has observed all of the regions forming high- mid-infrared flux threshold is used as a proxv for the pres- 

mass stars withm a distanc e of 3 kpc from the Sun, one of which ence/absence of a high-mass stellar embryo (see IMotteetalJ 

* Herschel is an ESA space observatory with science instruments pro- ^OOT), have short lifetimes. The high-mass class 0-like pro- 

vided by European-led Principal Investigator consortia and with impor- ^stars forming within IR-quiet MDCs ([Bontemps et al.| |2010D 

tant participation from NASA. See http://herschelesac.esa.int/ are observed to be fed on small scales by supersonic gas flows 

Figures 9-1 1 are only available in electronic form at ( Csengeri et al. 201 1). These short timescales and fast gas flows 

[http : //www, aanda. org are consistent with molecular clouds and dense fi laments dynam- 

^ \http:/ /hobys-herschel . cea.fr\ ically formed by colliding flows of H i gas (e.g. I Schneider et aO 
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l2Q10al: iNguyen Luong et al.ll201 ih . Cold and dense filamentary 




Fig. 1. The three-colour image built from Herschel images with 
red = 250 yum, green=160 yum, and blue=70 fim. The bright dif- 
fuse emission on the right corresponds to the Galactic plane. 
H II regions are prominent blue loops/bubbles, while earlier stage 
star-forming sites are red filaments. 

structures in molecular clouds are potential sites to find the pre- 
cursors of high-mass stars. One such example are Infrared dark 
clouds (IRDCs), which are dark extinction features against the 
Galactic background at mid-IR wavelengths. They have rather 
high column densities (>10^^ c m~^), col d temperatures (<20 K), 
and filamentary structures (e.g. iPeretto & Fuller 2010), bearing 
resemblances to the swept-u p gas features in t he simulations 
of co lliding atomic gas (e.g. .Heitsch et al .1120091: iBanerjee et al.l 
I2009h . 

The IRDC G035.39-00.33 (also called G035.49-00.30) is 
an IRDC filament lo cated in the W48 molecular complex at a 
distance of ~3 kpc (iRygl et al.l l2010h . Using ^^CO emission, 
ISimon et"aL 1 ^2006") estimate that the filament has a mass of 
-9000 Mq in an area with an eff^ectiv e radius of ~ 10 pc and a me- 
dian column density of 1x10^^ cm~^. |jimenez-Serra et al ] (12010) 
observed that extended SiO emission is associated with this fila- 
ment, which they interpreted as being produced by low- velocity 
shocks associated with colliding gas flows and/or shocks from 
protostellar outflows. Here we use Herschel data to investi- 
gate the star formation activity in G035. 39-00. 33 and determine 
whether any of the SiO emission could be associated with con- 
verging flows. 



2. Herscfie/ observations and ancillary data 

The entire W48 molecular complex was observed on 18 and 
19 September 2010 u sing PACS (iPoglitsch et al.1 2Q10) at 
70/160yum and SPIRE (iGriflan et al.l l2010h at 250/350/500yum 
in the parallel scan-map mode with a scanning speed of 20'7s. 
Two perpendicular scans were taken to cover a SPIRE/PACS 
common area of 2.5° x 2.5°. The data were reduced in two 
steps. The raw data (level-0) of the individual scans from both 



PACS and SPIRE were calibrated and deglitched using HIP^ 
pipeline version 7.0. The SPIRE and PACS level- 1 data were 
then fed to version 4 of the Scanamorphos software packageEl 
(Roussel 2011, submitted), which subtracts brightness drifts 
by exploiting the redundancy of observed points on the sky, 
masks remaining glitches, and produces maps. The final images 
have angular resolutions of ~6'', 12'', 18'', 25", and 37" and 
lo- rms of 0.02 Jy/ 1.4 "-pixel, 0.08 Jy/2.8"-pixel, 1 Jy/beam, 
1.1 Jy/beam, and 1.2 Jy/beam for 70yum, 160 yum, 250yum, 
350//m, and 500 yum, respectively. All Herschel images were 
converted to MJy/sr by multiplying the aforementioned maps 
by 21706, 5237, 115, 60, and 27 for 70yum, 160 yum, 250yum, 
350yum, and 500 yum, respectively. The final maps of W48 are 
shown as a three-colour image in Fig.[T]and individual images in 
Figs.|9t-e. 

To complement the spectral energy distributions (SEDs) 
of young stellar objects, we have used Spitzer data at 3.6, 
and 24 um (from the GLIMPSE and MIPSGAL surveys; 

" LABOCA data at 
Schuller et al] l2009lfl 



iBenjamin et al.l I2003L ICarey et al.1 1200911 
870 fim (from the ATLASGAL survey; 
and BOLOCAM data at 1.1 mm (from the BGPS survey; 
iBallv et al. 20100. We note that we used the images but not the 
compact source catalogues produced by these surveys since we 
extracted sources using a new algorithm simultaneously on all 
ten images (see Sect. 4.1). 

In the present study, we did not apply colour corrections for 
PACS, SPIRE, or ancillary data. These flux corrections would 
be rather small (<10% in the case of Herschel) and are covered 
by the 30% absolute calibration uncertainty we use when fitting 
SEDs (see Sect. 4.3). 



3. The G035.39-00.33 ridge characterized by 
Herschel 

W48 is a massive molecular cloud complex (8x10^ M© that has 
been identified in the extinction map of Bontemps et al. in prep.) 
located slightly south of the Galactic plane. The numerous H 
II regions (see Fig.[T]) indicate that this region has actively formed 
high-mass stars. Figure [T] also clearly displays future (high- 
mass) star-forming sites, such as the cold and dense filament 
IRDC G03 5.3 9-00.33, which is a prominent elongated struc- 
ture oriented north-south in the Herschel maps (see Figs.[2^-d). 
This site appears as a dark feature from Herschel 70 jim down to 
Spitzer 3.6 yum, and in emission from 160 jim upward. 

Prior to deriving any physical parameters, we added zero off'- 
sets, which were determined from the correlation with Planck 
and IRAS data following the procedure of Bernard et al. (2010) . 



^ HIPE is a joint development software by the Herschel Science 
Ground Segment Consortium, consisting of ESA, the NASA Herschel 
Science Center, and the HIFI, PACS, and SPIRE consortia. 

^ http://www2. iap.fr/users/roussel/herschel/ 

^ The Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE) and the Multiband Imaging Photometer for Spitzer 
GALactic plane survey (MIPSGAL) provide 3.6-24 fim images of the 
inner Galactic plane with 1.5'' - 18'' resolutions. Detailed information 
and reduced images are available at http://www.astro.wisc.edu/sirtf/ and 
http: //mips gal. ipac. caltech. edu/ 

^ The APEX Telescope Large Area Survey of the GALaxy 
(ATLASGAL) provides 870yL/m images of the inner Galactic plane 
with 19" resolution. Reduced images will soon be available at 
http: //www. mpifr-bonn. mpg. de/div/atlas gal/index, html 

6 The BOLOCAM Galactic Plane Survey (BGPS) provides 1.1 mm 
images of the inner Galactic plane with 33" resolution. Reduced images 
are available at http://irsa. ipac. caltech. edu/data/BOLOCAMjGPS/^ 
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Fig. 2. (a) PACS 70/im image of the G035.39-00.33 filament; (b) PACS 70yum image of the G035.39-00.33 filament and its 
surroundings; (c) SPIRE 250yum image of the G035. 39-00. 33 filament and its surroundings; (d) Temperature (colour) and column 
density (contours from 1.5 to 9 by 1.5 x 10^^ cm"^) images. The dense cores with mass >20 M© are indicated by black ellipses, 
those with mass <20 M© by white ellipses. The elliptical sizes represent the FWHM sizes at 160 fim. The extent of the IRDC 
(>3 X 10^^ cm"^) is indicated by a white polygon in b-d. 
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a(2000) a(2000) 

Fig. 3. Colum n density (a) and du st temperature (b) maps of W48 compiled from Herschel images with a HPBW of 37 '\ The dust 
opacity law of Hildebrand] (Il983h is used with a dust emissitivity (5 = 2. The W48 giant molecular cloud is indicated, the zoom on 
the IRDC G035. 39-00. 33 and its surroundings shown in Figs.[2]3-d is outlined. Plus signs (+) indicate the location of supernova 
remnants. 



We then derived the dust temperature and column density maps 
of W48 (Figs. [3^-b) by fitting pixel-by-pixel single grey-body 
SEDs. Only the four longer-wavelength Herschel bands, with a 
resolution of 37'' and equal weight, were used in our SED fit- 
ting, as the 70 yum data may not be tracing the cold dust that we 
are most interested in. We assumed that the emission is optically 



t hin at all wavelengths and took a dust opacity law of lHildebrandl 
(119831) with a dust emissivity index {3 = 2 (see equation given in 
Sect.O. 

The column density and temperature maps of Fig. [2li show 
that the IRDC harbours dense and cold material (A^H2 ~ 3 - 
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Table 1. Characteristics of the compact sources in G035. 39-00. 33 and its close surroundings: sources outside the IRDC are marked 
with an asterisk; FWHM deconvolved sizes are measured at 160yum; dust temperature and mass are derived from SED modelling 
(see Fig. 131); bolometric and submillimeter luminosities are measured from SEDs (see Sect. 14.3b and lower limit of ^>350/zm/^bol 
are given when sources are undetected at any wavelengths in the 3.6-70 micron range. 
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9 X 10^2 cm-2 and r^ust -13-16 K). With a length of ~6 pc 
and a width of -1.7 pc, this filamentary morphology is in good 
agreement with that found in the surface density map derived 
from 8 Jim extinction images (e.g. Butler & Tan 2009). The col- 
umn density of G035. 39-00. 33 derived from Herschel images is 
among the highest of known IRDCs ( Peretto & Fuller 2010) and 
is typical of "ridges", i.e. dominant filame nts with a high col- 
umn density as defined by lHill et"an (1201 ih . In contrast to most 
massive star-forming filaments, this IRDC has however a low 
temperature. From the column density map, we estimate IRDC 
G035.39-00.33 to have a total mass of -5000 M© within an area 
of ~8 pc^ corresponding to Ay > 30 mag as outlined in Fig.[2ll. 
This mass is similar to that derived from mid-infrared extinction 
dutler & Tan 2009) and agrees well with that measured from 
a column density map of th e IRDC where its background has 
been subtracted (method by iPeretto et al.|[2QTQl) . In the follow- 
ing, we show that such a cold ridge has a high potential to form 
(massive) stars. 

4. Massive dense cores in G035.39-00.33 and its 
surroundings 

4.1. Extracting Herschel compact sources 

Herschel compact sources were extracted using the multi-scale, 
multi- wavelength getsources algorithm (version 1.110614, see 



Im en'shchikov et al. and Men'shchikov in prep, for full de- 
tails). In our extraction, we first used only the five Herschel 
wavelength images for detecting sources and then all ten wave- 
lengths, from 3.6 jim to 1 . 1 mm, when measuring their fluxes. At 
the detection step, the five Herschel images were decomposed 
into multi-resolution cubes of single- scale single- wavelengt h 
images, as in the MRE-GCL method (see iMotte et al.ll2007h . 
The five cubes were then combined, with greater weight be- 
ing given to the higher resolution images, into a single cube of 
single-scale combined-wavelength images. The Herschel com- 
pact sources were finally detected within this cube, and their 
spatial positions and initial sizes were recorded. At the mea- 
surement step and in the initial ten wavelength images succes- 
sively, the sources properties (including FWHM size and inte- 
grated flux) were computed at their detected location, after the 
background had been subtracted and overlapping sources had 
been deblended. In the final getsources catalogue, each Herschel 
compact source has a single position and ten FWHM sizes at ten 
diff'erent wavelengths. 

We selected sources with deconvolved FWHM sizes at 
160 ^m smaller than 0.3 pc and a >7cr detection at more than 
two Herschel wavelengths (equivalent to >5cr measurements at 
more than six wavelengths), giving a catalogue of 28 sources 
complete down to 5 M© (see Figs. [2]3-c and Table [TJ. Among 
them, 13 sources were identified as probably the best high-mass 
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star precursor candidates (M > 20 M©), which we refer to 
as massive dense cores. We consider these 28 sources as robust 
sources because 92% of them were also identified by the cutex 
and csar algorithms ( Molinari et al.ll2Ql ll: Kirk in prep.). 

The sources studied here and more generally by the HOB YS 
key programme are cloud fragments that are strongly centrally 
concentrated but do not have firm outer boundaries when stud- 
ied from subparcsec to parcsec scales. Their inner part can be 
fitted by a Gaussian with a -0.1 pc FWHM size (e.g. lMotte et al.l 
|2007|), which can only be resolved by the PACS cameras at the 
0.7 - 3 kpc distances of HOBYS molecular cloud complexes. 
From submillimeter observations, these MDCs have been found 
to be quasi- spherical sources with a radial density distribution 
very close to the p(r) oc law for r ~ 0.1 - 1 pc (e.g. 



iBeuther et al .120021: Mueller et al.ll2002h . This density power law 
and thus a mass radial power-law of M(< r) oc r are theoreti- 
cally expected for dense cores formed by their own gravity, from 
singular isoth ermal spheres ( Shu 1977) to protocluster clumps 
(lAdamsll2000l) . For optically thin dust emission and a weak tem- 
perature gradient (T(r) oc r"^ with ^ = for starless or IR-quiet 
MDCs and q = 0.4 for IR-bright MDCs), one expects the fluxes 
integrated within the apertures to vary (close to) linearly with the 
angular radius: F(< 6) oc for IR-quiet protostellar dense cores 
and F(< 6) oc 6^-^ for IR-bright protostellar dense cores (see e.g. 
iMotte & Anddl200l| for more details in the case of protostellar 
envelopes). We thus expect that the deconvolved sizes and in- 
tegrated fluxes of the dense cores extracted here vary with the 
beam size and in practice increase from 160 to 500 yum, in line 
with the increase by a factor of >3 in the Herschel HPBW. This 
behaviour is indeed observed here (see Table 2), which we stress 
is independent of the compact source extraction technique used 
because it h as also been observed for sources extracted by the 
MRE-GCL (iMotte et al.ll2010 l) and the Cutex (Giannini et al. 
submitted) algorithms. 



Table 2. Scaled fluxes and quality of SED fits for the catalogue 
of Table. □ 



Wavelength [jim] 
HPBW 


160 
12 


250 
18 


350 
25 


500 
37 


< FWRMf" > [pcY 


0.18 


0.30 


0.40 


0.60 


^scaled ^ 
^original 




53% 


39% 


27% 


1 1 
l^original ~^original,GB | 

^original,GB 




90% 


220% 


510% 


\f f 

r scaled ^scaled,GB| 

^scaled,GB 




15% 


20% 


60% 



Notes: 

(a) : The deconvolved FWHM sizes at wavelength A were 
calculated from the average FWHM sizes measured by 
getsources and the beam sizes as 

FWRMf"" = V< FWHMx >^ -HPBW^ 

(b) : /^original ^scaled the original and scaled fluxes. 

(c) : i^original,GB is the grey-body flux fitted to the original fluxes. 

(d) : i^scaled,GB is the grey-body flux fitted to the scaled fluxes. 



4.2. Scaling Herschel fluxes for compiling SEDs 

To limit the influence of difl'erent Herschel resolutions and re- 
strict ourselves entirely to the MDCs/dense cores size (~Q-1 pc) , 
we followed the procedure introduced by iMotte et al.l (l2010l) . 



We kept the R\CS 70 and 160 jim fluxes unchanged and linearly 
scaled the SPIRE 250, 350, and 500 yum fluxes to the source size 
measured at 160yum. This process assumes that (1) the size mea- 
sured at 160yum reflects the spatial scale of the dense cores; (2) 
the 250, 350, and 500 yum emission are mainly optically thin; 
and (3) the F(<6) oc relation mentioned in Sect. 4.1 applies. 
In contrast to the >160yum emission, compact 70yum emission 
originates from hot dust close to the protostar and does not trace 
the cold component of a dense core. The relation F(<6) oc 6 
is correct if the density and temperature power laws mentioned 
above apply to a large portion of the dense cores, roughly from 
half the resolution at 160yum to the resolution at 500 yum, i.e. 
from 0.02-0.09 pc to 0. 1-0.6 pc depending on the cloud distance. 
However it does not preclude the subfragmentation, inner den- 
sity flattening, and/or inner heating of the (massive) dense cores. 
As explained in Sect. 4.1, these are reasonable assumptions for 
MDCs. In contrast, if an HII region develops, the density and 
temperature structures of its parent dense core will probably be 
strongly modified and such a scaling cannot apply. The Herschel 
fluxes (and other submillimeter measurements here) are scaled 
following the simple equation 
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Fax 



FWHM^ll 
FWHUf"" ' 



(1) 



where F^ and FWHUf"" are the measured flux and deconvolved 
FWHM size at wavelengths A > 250yum. For simplicity and 
because their characteristics do not impact our discussion, we 
kept a linear scaling for IR-bright MDCs #1, #2, and #3. 




MDC #17 

^13K 



J 



10 10 10 10 10 10 
Wavelength [ju-m] 

Fig. 4. SEDs built from Herschel and other wavelengths of: the 
UCH II region (a) and two MDCs associated with SiO emis- 
sion (b & c). The curves are grey-body models fitted to data 
at wavelengths > 160 yum. The single temperature grey-body fit 
(black curve) is consistent with the two temperatures grey -body 
fit (green curve). Error bars correspond to 30% of the integrated 
fluxes. 



The complete SEDs of the 28 dense cores of Table [T] are 
built from Herschel and ancillary data (see e.g. Fig.|4]). They can 
be compared to grey-body models if they arise from a unique 
gas mass. The scaling procedure used here allows for this, for 
a mass reservoir defined at 160yum. Grey-body models are not 
perfect models but closely represent the long wavelength com- 
ponent (160yum to 1.1 mm) of SEDs compiled for sources such 
as starless, IR-quiet, or IR-bright protostellar MDCs. For the 
example MDC #17, without scaling, Herschel provides almost 
constant fluxes from 250yum to 500 yum (see Fig.O, making the 
grey-body fit very unreliable. A few other SEDs even increase 
from 160 yum to 500 yum, unrealistically suggesting a <5-10 K 
dust temperature. Without complementary submillimeter obser- 
vations of high angular resolution and without any scaling, the 
SED of several cloud fragments, such as MDC #17, could not 
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Fig. 5. Example SED compiled from fluxes with (in black) and 
without (in red) scaling for MDC #17. The SED with scaled 
fluxes portrays a more classical behaviour for the Ray leigh- Jeans 
part of the protostellar dense core. 

be fit at all. As shown in Fig. \5\ the scaled SPIRE fluxes por- 
tray a more classical behaviour for the Rayleigh- Jeans part of 
a protostellar dense core SED. As can be seen in Table [2l the 
250 Jim, 350 yum, and 500 jim fluxes are quantitatively lower but 
these flux changes help to improve the SED fit by a grey-body 
model. Indeed, the variations between the measured Herschel 
fluxes and the fitted grey-body fluxes are -5-10 times larger in 
the case of the original fluxes than for scaled fluxes (see Fig.O. 
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Fig. 6. The absolute value of the normalized difl'erence between 
the original fluxes with their best-fit grey-body fluxes (crosses) 
and between the scaled fluxes with their best-fit grey-body fluxes 
(circles). 

After flux scaling, the dust temperatures obtained from the 
grey-body models of sources within G035.39-00.33 are dis- 
tributed around <T^^^{>^ 14.5 K (see Fig. [7^). They are in closer 
agreement with those measured from the dust temperature im- 
age (-16.5 K, see Fig. [2li and [7^) and are -20% higher than 
those obtained from grey-body models for the original fluxes 
(-11 K). Most dust temperatures derived from SED fits to the 
original SPIRE fluxes are too cold. We do expect the temper- 
ature of compact cloud fragments to difl'er from that of their 
surrounding cloud and in practice to be lower if the protostel- 
lar dense cores are barely evolved as in IR-quiet massive dense 
cores. However, a -5 K drop in temperature is quite extreme 
when averaged over only a three times smaller cloud fragment 



and a <T^^^{>-- 11.3 K for MDCs is unusually cold. For in- 
stance, a temperature of -20 K has been measured in Cygnus X 
massiv e dense cores from ammonia measurements and Herschel 
SEDs dMotte et al.ll2007l: Hennemann et al. in prep.). 
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Fig. 7. Dust temperature (a) and mass (b) distributions of the 
28 compact sources derived from fitting a grey-body model to 
the original (in black) and scaled (in red) fluxes at wavelengths 
/I > 160 yum. The green histogram is the temperature distribution 
of the filament IRDC GO.35.39-00.33 taken from Fig.[2]l. 



4.3. Deriving physical parameters for the compact sources 

After scaling the fluxes, we compiled a SED and fit a single 
grey-body to the 160yum-l mm wavelength data (see Figs. S^- 
c and [TTh-e) to derive the mass and the average dust tem- 
perature of each dense core. We assumed a dust opacity law 

kq X (30^^) ^ with kq = 0.1 gcm^ and a dust emissiv- 



ity index ofjS = 2. Our limited knowledge of the dust emissivity 
makes masses systematically uncertain by a factor of -2. The 
mass and dust temperature uncertainties quoted in Table [T]reflect 
those of SED fits to fluxes varying by 30% with a constant emis- 
sivity. We successfully counterchecked our results obtained with 
single grey-body fits to the long wavelength (>160yum) compo- 
nent of SEDs by fitting two-temperature grey-body models to 
the data (see Fig. [4]). 

We also derived the bolometric luminosity and submil- 
limeter luminosity L^ssojum by integrating fluxes below the mea- 
surement points i.e., from 3.6 jim to 1 . 1 mm and from 350 jim to 
1.1 mm, respectively. When sources were not detected at mid- 
infrared wavelengths. Table [T] gives a range of bolometric lumi- 
nosities using the 3.6, 8, and 24 jim noise levels as upper limits 
or setting them to 0. 

Cloud fragments in our sample have masses ranging from 
1 M© to -50 M© with deconvolved FWHM sizes at 160 jim rang- 
ing from 0.1 to 0.3 pc and temperatures from 10 to 30 K. The 
mass of the gas reservoir selected at 160yum and measured by 
a grey-body fit is, on average, 40% that obtained from the SED 
fit of original fluxes (see Fig. (TJ)). The shrinking, by a factor of 
-2, of the mass reservoir arises from the decision to focus on the 
gas contained within the 160 jim sizes, which are twice as small 
as the mean sizes at the SED peak (generally close to 250yum). 
Without such a scaling, the present study would have considered 
cloud fragments that are, on average, twice as massive but also 
-8 times less dense. These fragments would thus not correspond 
as well to the typical scale and density of MDCs and would logi- 
cally have a less efficient transfer of matter from the clump gas to 
star. For the reasons given in Sects. 4.1-4.2 and despite the modi- 
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fications implied, we estimate that the rescahng of SPIRE fluxes 
is a necessary step in a simple SED fitting approach applied to 
large samples of compact sources. 



5. Discussion 

5.1. A burst of star formation activity 

We have focused here on the 13 best high-mass star pre- 
cursor candidates, or MDCs, which are the most massive 
(M > 20 Mq) and dense (FWHM averaged density of 
>2 X 10^ cm"^) cloud fragments, and thus the most compact 
(FWH M < 0.19 pc). Several M DCs lie within the clumps identi- 
fied by lRathborne et al.l (l2006l) at 1.3 mm, which are larger (0.2- 
0.9 pc) cloud structures. 

The MDC #1 was classified as an Ultra-Compact H ii 
(UCH ii) region bas ed on its free-free emission at 1 .4 and 5 GHz 
( Becker etal.lll994 . The MDC #2 is an UCH ii candidate ac- 
cording to the Red MSX Source survey (iMottram et al.l l201 ih 
but no centimeter free-free emission has yet been dete cted. Water 
maser emission was detected towards MDC #2 ( Cha mbers et al.l 
120091) . These two MDCs together with MDC #3 are the most 
luminous MDCs in our sample, i.e, L^qi ~ 3000 - 5000 L©, 
while others have L^qI < 300 L©. The high luminosity of MDCs 
#1, #2, and #3, their warm dust and high mass imply that they 
are already evolved massive star-forming sites such as IR-bright 
MDCs (#2 and #3) or UCH iis (#1). These three MDCs also 
have distinct 24yL/m fluxes (>1 Jy) compared to the remaining 
cores (<1 Jy). This flux limit is close to that used to define IR- 
quiet MDCs (see |Motte_etaD '2007 after distance correction). 
The separation based on 24 fim fluxes is also consistent with the 
L>35o^m/^bol ratio, which is small (< 1%) for IR-bright MDCs 
or UCH II regions and large (> 1-20%) for IR-quiet MDCs har- 
bouri ng high- mass class Os or for starless dense core candi- 
dates. lAndre et al.l (|2P00) used this ratio to separate low-mass 
class and class I protostars (size ~ 0.01 - 0.1 pc) with 1% 
as the dividing line. Since MDCs are larger-scale cloud struc- 
tures (~ 0.1 pc), the ^>350//m/^bol ratio of IR-quiet MDCs is ex- 
pected to be substantially larger than those of low-mass class Os. 
Among the remaining ten MDCs with low 24 fim emission, one 
(MDC #26, outside the G035.39-00.33) is a good candidate star- 
less dense core since it has no emission at wavelengths shorter 
than 70 jim. 

Interestingly, 55% (15) of the dense cores and 70% (9) of 
the MDCs are located inside the filament G035.39-00.33, which 
has a mass density of 600 M© pc"^ and a mass per unit length 
of 800 Mq pc"^. It is much higher than those of the highly ex- 
tincted regions ( Ay > 30 mag) of Perseus and Ophiuchus (see 
iHeiderman et al.t .2010). This high concentration of (massive) 
dense cores indicates that the gravitational potential of G035.39- 
00.33 helped them to build up. From their SED, the MDCs 
in the G035. 39-00. 33 filament are IR-quiet protostellar dense 
cores and thus harbour young protostars (see Table [T]). This 
predominance of IR-quiet MDCs suggests that most of these 
dense cores have been formed simultaneously and are probably 
forming just following the formation of the filament. In agree- 
ment wit h the presence of an UCH ii region and previous stud- 
ies iMotte et al.1 120071: lRusseil et al.1 120101) . we assumed a 
20 - 40% mass transfer from MDC to stellar mass. We there- 
fore estimate that between four and nine of the MDCs (those 
with >20 - 40 M©, see Table 1) would indeed form high-mass 
stars in G035. 39-00. 33. With nine high-mass protostars form- 
ing in this 5000 M© fi lament and assuming the initial mass func- 
tion of lKroupal (l200ll) . the star formation eflftciency (SEE) could 



be as high as -15%. The numerous protostars (being detected 
at 70yL/m) observed within the G035. 39-00. 33 filament (area of 
~8 pc^) shown in Eig. [2^, are estimated to have a mean mass 
on the main sequence of ~2 M© (meaning that stars form in 
5-50 Mq dense cores with 25% mass transfer efficiency) and a 
lifetime of ~10^ yr, close to the dense-core free-fall time. This 
implies a star formation rate (SER) of ~300 M© Myr"\ which 
exc eeds the mean SER of Gould Belt clouds by a factor of seven 
(see lHeiderman et al. 2010) a nd is close to that of the W43 mini - 
starburst region (MoUq et al1l2003l: iNguven Luong et al.|[2QTTh . 
The SER density within this ~8 pc^ high column-density (Ay > 
30 mag) filament, I^sfr ~ 40 M©yr"^kpc"^, is high and has 
only been measured in Gould Belt dens e clumps (10-100 times 
smaller areas with Ay > 20-40 mag, see IHeiderman et ani2010l: 
iMaury et"aD 20TT) and in ex tragalactic circumnuclear starburst 
regions (e.g.^ennicutt 19981) . 

The measurements (SEE, SER, I^sfr) of this high star for- 
mation activity are uncertain by a factor >2, coming essentially 
from the uncertainty in the number of high-mass protostars and 
the mass transfer efficiency from dense cores to stars (see above). 
In contrast, the flux scaling should have a relatively low influ- 
ence on our results. Indeed, without scaling, massive clumps are 
twice as massive but would form stars of the same stellar mass 
on the main sequence. Despite these limitations, the SEE, SER, 
and Us FR values above suggest that a mini-starburst event, i.e. a 
miniature model of events in starburst galaxies, is occurring in 
this filament. To keep the SEE in the entire W48 complex con- 
sisten t with those seen in other clouds (1-3% per 10^ yrs, ISilkl 
Il997h and reconcile the SEE and SER values given above, the 
star formation event should be short (~10^ yr). If confirmed, this 
mini-burst of recent star formation would be consistent with the 
filament G035.39-00.33 being formed through a rapid process 
such as converging flows. Another interpretation could be that 
star formation has recently been triggered, but a nearby trigger- 
ing source has not been found. 



5.2. Are protostars responsible for the SiO emission? 



iJimenez-Serra et al.l ( |2010l) mapped the filament G035.39-00.33 
in SiO (2-1) (HPBW^ 28'', see Eig. [8]) finding 0.01 - 
0.05 Kkms"^ widespread emission along part of this filament. 
They proposed that it could be produced by low- velocity shocks 
associated with converging flows and/or outflows from proto- 
stars. SiO is indeed an excellent tr acer of shocked gas associ- 
ated with molecular outflows (e.g. ISchilke et al.l 1 19971) and is 
particula rly prevalent towards early-stage hi gh-mass protostars 
(e.g. Hat cheU et ani200ll: iMotte et al.l l2007l) . It has also been 
proposed that low-velocity shocks in a high-density medium 
could sputter dust grains and enhance SiO in the gas phase 
(Lesafl're, R, Gusdorf, A. priv. comm.). 

We investigated the origin of this SiO emission by correlat- 
ing its spatial distribut ion w ith that of young protostars detected 
by Herschel (see Sect. 15. II and Table [T]). Two of the three peaks 
of SiO (North and South) coincide with IR-quiet massive dense 
cores: MDCs #6 and #17 (see Eigs. [8]and[T0|). Given that the 
SED of MDC #17 displays a dominating cold and massive enve- 
lope and a warm component detected with Spitzer (see Eig. lit), 
it most probably hosts an early-stage high-mass protostar. The 
MDC #6 is not detected at wavelengths shorter than 24 [im (see 
Eig. lit) but its compact 70 yum emission is three times that of the 
flux predicted by the grey-body model that fits its cold envelope. 
This MDC also c oincides with a 40'' -beam water maser source 
(Ry gl etal.ll2QTl) suggesting that it harbours an intermediate- 
or high-mass protostar. The SiO emission at these two positions 
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Fig. 8. (a) Part of G035.39-00.33 seen at 70 jim (colour) and in 
SiO (contours from 0.01 to 0.05 K-km s"^ by 0.01 Kkm s"^ from 
iJimenez-Serra et al]|2010h . The dense cores with mass >20 M© 
are indicated by red elHpses, those with mass <20 M© by 
white elHpses. (b) Zoom towards source #27 seen at 8-160 jim. 
Red and black boxes indicate the 8 yum-only and 24 yum-only 
sources, respectively. 

could easily originate from shocks within protostellar outflows 
because they are comparable to the weakest ones a ssociated with 
IR-quiet MDCs in Cygnus X (iMotte et al.ll2007l after distance 
correction). 

In contrast, the eastern SiO peak, although detected in 
SiO (2-1) at -0.04 Kkm s"\ does not coincide with any reliable 
protostar candidate. Source #27, which lies within the SiO beam, 
is only detected shortwards of 160 yum (see Fig.lHJ)). A specific 
run of getsources, allowing low-reliability measurements (a 3(T 
detection at a single Herschel wavelength), has in fact been done 
at the location of the eastern SiO peak to extract source #27. 
The SED of this source (see Fig. \TTti ) consists of two unrelated 
components: one mid-IR source of unknown nature detected 
at 8, 24, and 70 yum plus some nearby cloud structure whose 
160 yum-1.1 mm emission extends until the location of source 
#27. Re-examining the Herschel and Spitzer images, we estimate 
that source #27 could be at most a 1 M© dense core containing 
an evolved low-mass class I protostar. Since SiO outflows are 
known to be rare and confined to the proximity of interm ediate- 
mass class and high-niass protostars (e.g. Martin-Pintado et al.l 
119971: iGueth et al.lll998l: iMotte et al.ll2007h . there is only a low 
probability that such a protostar could drive a strong-enough 
outflow to be responsible for the SiO emission of the eastern 
peak. Two other sources observed within the SiO beam are only 
detected at 8yum or 24/im and are thus probably even more 



evolved young stellar objects or variations in the background 
emission. Outflows from neighbouring low- to high-mass pro- 
tostars, >0.7 pc away from the eastern SiO peak, would also 
be unable to create such a SiO emission. Indeed, the extreme 
SiO outflow arising from the intermediate-mass star LI 157 has 
e mission peaks th at are at less than 0.2 pc from the protostar 
(Gueth e t aDll998l) . The protostars here are unlikely to be ex- 
treme and less confined. 

Since no definite protostar is detected by Herschel towards 
the eastern SiO peak, the interpretation that the SiO emission 
could partly originate from large-scale converging-flows may 
be valid, at least at this location. Higher-resolution and higher- 
sensitivity observations would be necessary to definitively reach 
this conclusion. Further studies of the SiO emission and its rela- 
tion to other high-density gas tracers, as well as the dust proper- 
ties, will be presented by Henshaw et al. (in prep). 

6. Conclusion 

We have used PACS and SPIRE maps of Herschel to investigate 
the star formation activity in the IRDC G035.39-00.33, a cold 
(13 - 16 K) and dense (A^H2 up to 9 x 10^^ cm"^) filament, which 
we have qualified as a "ridge". We have proposed a new ap- 
proach to analysing the SED of compact sources compiled from 
Herschel fluxes. We fitted single grey-body models to fluxes 
that had been linearly scaled downwards to correspond more 
closely to the mass reservoir emitting at 160 yum, the wavelength 
at which emission from the cold envelope of compact source is 
probed with highest angular resolution. On the basis of this pro- 
cedure, we have found that Herschel detected a total of 28 dense 
cores (FWHM^O.15 pc), among them 13 MDCs with masses of 
20-50 M© and densities of 2-20 xlO^ cm"^, which are poten- 
tially forming high-mass stars. Given their concentration in the 
IRDC G03 5.3 9-00.33, they may be participating in a mini-burst 
of star formation activity with SFE-^ 15%, SFR-^ 300 M© Myv-\ 
and lis FR ~ 40 M© yr"^ kpc"^. Two IR-quiet MDCs could be the 
origin of most of the extended SiO emission observed, the re- 
mainder possibly originating from a low- velocity shock within 
converging flows. 
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Q. Nguyen Lifo'ng et aL: Massive star formation in G035. 39-00.33, a mini-starburst event 




a(2000) a(2000) 




a(2000) 

Fig. 9. Herschel images of W48: a) PACS 70 jim (HPBW ~ 6''), b) PACS 160 jjm (HPBW ~ 12'0, c) SPIRE 250 jim (HPBW ~ 
18'0, d) SPIRE 350 fim (HPBW ~ 25''), and e) SPIRE 500 jim (HPBW 37''). The bright diffuse emission on the right of each panel 
is from the Galactic plane. Plus signs (+) indicate the location of supernova remnants. 
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Fig. 10. Zooms towards MDC #6 (upper rows) and MDC #17 (lower rows) at 8 - 250 fim. Cloud fragments with mass <20 M© are 
mar ked as white and MDCs as black ellipses. SiO contours are the same as in Fig. [8] (from 0.01 to 0.05 Kkms"^ by 0.01 Kkms"^ 
from I Jimenez- Serra et al 
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Fig. 11. SEDs built from Herschel and other wavelengths for sources lying towards the SiO peaks (except for those shown in 
Fig.lHl. The curves are grey-body models fitted for data at wavelengths > 160 jim. The one-temperature grey-body fit (black curve) 
is consistent with two-temperature grey-body fits (green curve). Error bars correspond to 30% of the integrated fluxes. No model 
can reasonably fit source #27. 



